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We show that, up to multiplication by constants, a Gaussian 
process has an infinitely divisible square if and only if its covariance 
is the Green function of a transient Markov process. 

1. Introduction. The question of the infinite divisibility of squared Gaus- 
sian vectors is an old problem which was first raised by Paul Levy in 1948 
[10]. Given a centered Gaussian vector (^>i, . . . , <p p ), when can the vector 
((j>i , . . . , (ftp) be written as a sum of n i.i.d. p vectors for every n £ N? Many 
authors have worked on this problem. We refer the readers to [6, 8, 13, 14, 17] 
and the references therein for more on this problem. In 1984, Griffiths [9] 
established a characterization of the p-dimensional centered Gaussian vec- 
tors with an infinitely divisible square. His criterion is difficult to use since 
it requires the computation of the signs of all the cofactors of the covari- 
ance matrix. Indeed, except for the Brownian motion (and, more generally, 
Gaussian Markov processes), there were no examples, in the literature, of 
processes satisfying this remarkable property nor examples of processes lack- 
ing it. 

We have recently shown [7] that the family of fractional Brownian motions 
provide examples of both kinds. A fractional Brownian motion is a real- 
valued centered Gaussian process with a covariance given by 

g(x,y) = \xf + \yf - \x-yf, 

where the index (5 is in (0,2). We proved that when (3 is in (0,1], then 
the square of this process is infinitely divisible, and when (3 is in (1,2), it 
is not. The critical index 1 corresponds to the Brownian motion. In both 
cases we have used the criterion of Griffiths. We have shown in [7] that if 
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the Green function of a transient Markov process is symmetric, then it is 
the covariance of a Gaussian process with an infinitely divisible square. In 
particular, when the index (3 is in (0, 1] , the covariance of the corresponding 
fractional Brownian motion can be interpreted as the Green function of the 
symmetric stable process with index ((3 + 1), killed at its hitting time of 
0. To treat the second case, we have shown directly that the condition of 
Griffiths is not satisfied. 

In view of these examples, a natural question arises. Is the representation 
of the covariance function of a centered Gaussian process, as the Green func- 
tion of a symmetric Markov process, a necessary condition for the infinite 
divisibility of its square? We show in this paper that, up to a multiplication 
by a constant function, the answer is affirmative. The result is presented in 
Section 3 in the form of a necessary and sufficient condition. The proof is 
based this time on a criterion for infinite divisibility established by Bapat 
[1], which we recall in Section 2. Although, at first sight, the verification of 
this criterion is as difficult as that of the equivalent criterion of Griffiths (one 
has to check here, too, the sign of each cofactor of the covariance matrix), 
it allows to significantly shorten the arguments. 

Gaussian processes with a covariance equal to a Green function play an 
important role in the study of Markov processes. Indeed, the Isomorphism 
theorem of Dynkin [5] provides an indentity in law connecting each of these 
Gaussian processes to the local time process of the corresponding symmet- 
ric Markov process. This identity has been exploited to study properties of 
the local time process of the Markov process, using similar properties of 
the Gaussian processes. We choose to mention here only two references, but 
many more papers on the subject are quoted in [2], for example. In [11] Mar- 
cus and Rosen have studied sample path properties of the local time process 
using similar properties of the Gaussian process. In [2] Bass, Eisenbaum and 
Shi used the Isomorphism theorem to establish the transience of the most 
visited sites of a symmetric stable process. The question of characterizing 
the Gaussian processes with a covariance that is equal to a Green function of 
a symmetric Markov process has been open since the Isomorphism theorem 
of Dynkin was first proved. The importance of an answer to this question 
is twofold; it will give a powerful tool for the study of these Gaussian pro- 
cesses, as well as that of their associated Markov processes. The results of 
Section 3 provide an answer to this question. 

In Section 5 we show that the Brownian sheet does not have an infinitely 
divisible square. Recalling that linear combinations of the squared com- 
ponents of centered Gaussian vectors are infinitely divisible, we find this 
example and that of the fractional Brownian motion, with index in (1,2), 
somewhat counterintuitive. 
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2. The criterion of Bapat. Let A = (^4ij)i<ij< p be a p x p matrix. We 
write A > if A tj > for all i,j. 

Definition 2.1. A matrix A = (Aij)i<ij< p is said to be an M-matrix 
if we have the following: 

(i) Aij < for all i^j, 

(ii) A is nonsingular and A^ 1 > 0. 

We refer the reader to [3] for the theory of M-matrices. 

Definition 2.2. A diagonal matrix is called a signature matrix if all its 
entries are either 1 or —1. 

The Laplace transform, ip(ti, . . . ,t p ) of the square of a p-dimensional 
Gaussian vector is given by 

V>(t) = [det(/ + GT)]- 1 / 2 , 

where t = (t±, . . . ,t p ), G is a positive definite p x p-matrix, T is the diagonal 
matrix with entries T„ = ti , and / is the p x p-identity matrix. 

Bapat [1] obtained the following characterization of the matrices G for 
which the Laplace transform ifi is infinitely divisible, that is, for which t/) S (t) 
is a Laplace transform for any <5 > 0. 

Theorem A. The Laplace transform ip is infinitely divisible if, and only 
if, there exists a signature matrix S such that SG~ 1 S is an M-matrix. 

Remark 2.3. It is elementary to check, using Theorem A, that a cen- 
tered three-dimensional Gaussian vector (771,772;%) such that 

E(7 ?1 7 ?2 )E(7 ?2 7 ?3 )E(7 ? 37 ? l)<0 

cannot have an infinitely divisible square. 

3. A necessary and sufficient condition for infinite divisibility. In a pre- 
vious work [7] on squared Gaussian processes, we have established the fol- 
lowing result. 

Theorem B. Let g be the Green function of a strongly symmetric tran- 
sient Borel right Markov process with state space (E,£). Let n be a centered 
Gaussian process with covariance g. Then the process rj 2 is infinitely divisi- 
ble. 
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The set E (by the assumption of a right process) is a Borel subset of 
a compact metric space and £ is the cr-field of its Borel sets. Note that if 
(r/ 2 (x),x 6 E) is an infinitely divisible squared Gaussian process, then, for 
any ^-measurable real valued function d, the process (d 2 (x)r] 2 (x), x E E) is 
also infinitely divisible. We therefore have the following: 

Corollary 3.1. Let g be the Green function of a symmetric transient 
Markov process on a state space (E,£). Then for any £ -measurable real val- 
ued function d on E, there exists a centered Gaussian process r\ with covari- 
ance equal to (d(x)g(x,y)d(y),(x,y) £ E x E). The process rj 2 is infinitely 
divisible. 

We shall first treat the case when E is a finite set. Theorem 3.2 completes 
Corollary 3.1, for this case, by showing that its sufficient condition on the 
covariance is also necessary for the infinite divisibility of rj 2 . 

Theorem 3.2. Let n be a centered Gaussian vector, indexed by a finite 
set E, with a positive definite covariance function (G(x,y), (x,y) 6 E x E). 
The vector rj 2 is infinitely divisible if, and only if, there exists a real valued 
function d on E such that, for any x,y € E, 

G(x,y) = d(x)g{x,y) d(y), 

where the function g is the Green function of a transient symmetric Markov 
process. 

Remark 3.3. Note that the Green function g of a symmetric Markov 
process X is always positive definite. Indeed, it is semi-positive definite 
(see, e.g., [11] or [7]) and it has been shown in [7], Section II, that, for 
any x±,X2, ■ ■ ■ ,x n in the state space of X, the matrix (g(xi,Xj))i<ij< n is 
invertible. 

Let r\ be a centered Gaussian process indexed by an infinite set E. Then 
it has an infinitely divisible square if, for every finite subset F of E, the co- 
variance of (ri x ,x 6 F) satisfies the condition of Theorem 3.2. This does not 
guarantee, a priori, that, for some deterministic function d, the covariance of 
(d(x)r] x ,x € E) is the Green function of a transient symmetric Markov pro- 
cess. Restricting our attention to E = M, and under the additional continuity 
assumption on the covariance function, the following theorem establishes the 
necessity of that condition. 

Theorem 3.4. Let n be a centered Gaussian process, indexed byM., with 
a positive definite covariance function (G(x,y), (x,y) £ M. 2 ). Assume that G 
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is jointly continuous. If the process rj 1 is infinitely divisible, then there exists 
a measurable real valued positive function d onR such that, for any x, y £ R, 

G(x,y) = d(x)g(x,y) d(y), 

where the function g is the Green function of a symmetric transient Markov 
process. 

Theorem 3.4 is proved in Section 4. Although Theorem 3.2 was a good 
hint to anticipate Theorem 3.4, we could not use it directly to prove it. 
Our argument is based on an explicit construction of the function d and 
of the Green function g of the Markov process that will be associated with 
the Gaussian process n. As a first step of this construction, we show that 
the covariance G has to be positive. The proof of Theorem 3.4 remains 
valid, with simple changes, when the index set R is replaced by a separable 
locally compact metric space. Moreover, we will see in Section 4 that, as a 
by-product of the proofs of Theorems 3.2 and 3.4, we obtain the following 
characterization of the associated Gaussian processes, which makes use of 
the definition below. 

Definition 3.5. A p x p matrix A is said to be diagonally dominant if, 
for every i = l,...,p, Yij=i A ij > °- 

Theorem 3.6. (i) A positive definite matrix G is the Green function 
of a finite state space transient Markov process, if, and only if G~ l is a 
diagonally dominant M -matrix. 

(ii) Let E be a separable locally compact metric space. Let (n x ,x £ E) be 
a Gaussian process with a continuous positive definite covariance (G(x,y), 
(x,y) G E 2 ). Then G is the Green function of a transient Markov process 
on E, if and only if, for every x±,X2, ■ ■ ■ ,x p in E, the inverse of the matrix 
(G(xi,Xj))i<ij< p is a diagonally dominant M-matrix. 



4. Proofs. 



Proof of Theorem 3.2. Let E be the finite set {xi, . . .,x p }. Theorem 
A guarantees the existence of a signature matrix S with diagonal Si,i = 
l,...,p such that the covariance matrix of (SiV Xi ,i = l,...,p) is positive. 
We shall therefore assume, from the onset, that the covariance G is positive. 

We will actually prove that the covariance G satisfies 

(1) G(x,y) = ^g(x,y), 

where the function g is the Green function of a transient Markov process 
and D is a strictly positive function. 
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This will be sufficient to establish Theorem 3.2. Indeed, if G satisfies (1), 
then we have 

g(x,y) = ^7~y G ( x >^)- 

Let U be the corresponding potential with density g(x,y) with respect to a 
reference measure [i. That is, 

Uf(x) = J g(x,y)f(y)n(dy). 

We then set m(dy) = n(dy)/D 2 (y). With respect to m, the kernel U has 
densities g(x,y) = D(x)G(x,y)D(y). Consequently, g is a symmetric Green 
function. Setting d = 1/D, one obtains the necessity of the condition of 
Theorem 3.2. 

To prove (1), let G be the covariance matrix of a p-dimensional centered 
Gaussian vector with an infinitely divisible square. By Theorem A, G^ 1 is 
an M-matrix. This implies (see [3]) that 

(2) G~ 1 = cl-B, 

where B > and c is strictly greater than the absolute value of any eigen- 
value of B. Further, by [3], Chapter 6, page 137 M36, since G" 1 is an M- 
matrix, there exists a diagonal matrix D such that Da > for all i and 
DG D has strictly positive row sum. This means that, for any i, 

According to Definition 3.5, DG~ 1 D~ 1 is diagonally dominant. Set T = 
D(\B)D- 1 . Then 

DG^D- 1 = D(cl - B)D~ 1 = c(I - T). 
Note that, for any i, 

The matrix T is therefore the transition matrix of a transient Markov chain 
(X n ) ng N with state space E = {x\ , X2, ■ ■ ■ , x p } satisfying Tij = F Xi {X\ =Xj), 
and YJj=i (-Xi = Xj) = 1 — ¥ Xi {X\ = A) , where A denotes a cemetery 
point. Let denote the total number of visits of X to the state xj. The 
Green function of X is defined by 



g(xi,Xj) =E a 
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It can be computed as follows: 

/ OO \ oo 

E Xi (&) = E Xi £ i {Xn=a; . } = £ = ( J " T)r.i, 

\n=0 / n=0 

which is defined and is finite, since, by (2) and the discussion following it, 
the spectral radius of T is strictly smaller than 1. Hence, we can write 

cDGD~ l = g; 

that is, for every x, y in E, 

cG(x,y) = ^g(x,y). 
D(x) 

We shall now use the well-known method to turn a finite state space Markov 
chain into a continuous time Markov process by subordination to a Poisson 
process with rate c (see, e.g., [4]). By its construction, this Markov process 
is transient with potential density (Green function) equal to g and (1) is 
established. □ 



Remark 4.1. Let Y be the Markov process with potential density g. 
Equation (1) looks as if G is the potential density of an h-path transform 
of Y . This is really the case if D is excessive for Y. Unfortunately, this is 
not true in general. Consequently, we see that the collection of covariance 
functions that correspond to Gaussian processes with an infinitely divisible 
square is somewhat richer than the set of Green functions of symmetric 
Markov processes. This remark remains true also when the index set E of 
the Gaussian process is infinite. 

To select covariance matrices that correspond to Green functions of Markov 
processes, we have the condition given by Theorem 3.6(i). Indeed, assume 
that G _1 is a diagonally dominant M-matrix. Keeping the notation of the 
proof of Theorem 3.2, we can choose D = I and obtain T = -B. For any i, 
J2*j=i Tij < 1- Since the spectral radius of T is strictly smaller then 1, for at 
least one i, the above inequality must be strict. Therefore, T is the transi- 
tion matrix of a transient Markov chain and we can conclude that G is the 
Green function of a transient Markov process with finite state space. 

To see that the condition is also necessary, consider the Green function 
g of a transient symmetric Markov process X with a finite state space E = 
{x\, X2, ■ ■ ■ , x p }. Then the inverse of the matrix G= (g(xi,Xj))i<i t j< p is a 
diagonally dominant M-matrix. Indeed Gij = Xj(I — T)^ , where Xj is the 
expected value of the exponential sojourn times in state j, T is the transition 
matrix of the Markov chain (X(5 n )) ne N) and S n is the nth jump time of X. 
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Proof of Theorem 3.4. We first show that G has to be positive. For 
a fixed n, we define the function d n on R by 

t/\ ^ • r k ^ fe "f" 1 

d n [x) = — it — < x < . 

Let K be the compact set [a, b] with a <b, and iet d n (K) = {d n (x) : cc € K}. 
The set d n {K) is finite set. Since the process (Vdn(x)> x e -^0 * s infinitely 
divisible, thanks to Theorem A, there exists a signature function s n (i.e., a 
function taking values in { — 1,+1}) on d n {K) such that, for every x,y € K, 

s n (d n (x))s n (d n (y))G(d n (x),d n (y)) > 0. 

Note that 

(3) s n (d n (x))s n {d n {y))G{d n (x),d n {y)) = \G(d n (x),d n {y))\. 

Since G is continuous, we obtain, by letting n tend to oo, that ]im n (s n (d n (x)) x 
s n (d n (y))) exists for all (x,y) for which G(x,y) / and, for such (x,y), 

(lim s n (d n (x))s n (d n (y))\G(x,y) = \G(x,y)\. 

For (x,?/) in K 2 such that G(x,y) = 0, there exists a finite sequence ax,a,2, 
...,a p such that G(x,ai)G(a±,a2) • • •G(a p -i,ap)G(a p ,y) ^ 0. Indeed, for 
2; 6 if, set C(z) = {u S R : u) ^ 0}. For each z, C(z) is an open set and 
\J ze K^( z ) i s a covering of the compact set K. Thus, there exists a finite sub- 
covering C(zi),C(z2), ■ ■ ■ ,C(z m ) of the set K. The sets of the covering are 
not disjoint. Let C(zi 1 ) be one of these sets. C(zi 1 ) is a union of disjoint inter- 
vals. If C(zi 1 ) does not cover K, then there exists Zi 2 in {zi,Z2, ■ ■ ■ , z m } such 
that C(zi 2 ) covers some of the endpoints of C(zi 1 ) that are in K. If C(zi 1 ) U 
C(zi 2 ) do not cover K, there exists Zi 3 in {z±,Z2, ■ ■ ■ ,z m } such that C(zi 3 ) 
covers some of the endpoints of C(zi 1 ) U C(zj 2 ) that are in K, and we may 
continue on until we exhaust all the finite covering above. Then we just have 
to make use of the sequence (z^ , Zi 2 , . . . , Zi m ) to construct (a\, ct2, . . . , a p ) (p< 
m), connecting x to y such that G(x, ai)G(ai, 02) • • • G(a p _i, a p )G(a p , y) ^ 0. 
Since 

Sn{d n {x))s 2 n {d n {ai))s 2 n (d n (a2)) ■ ■ ■ s 2 n {d n {a p ))s n (d n {y)) = s n (d n {x))s n {d n {y)), 

we obtain, using (3), the existence of lim n s n (d n (x))s n (d n (y)) for all x, y. 
Set 

H(x,y) = lim s n (d n (x))s n (d n (y)). 

The function H is symmetric and, by (3), 

H(x,y) = siga(G(x,y)), 
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and by its definition for all x,y,z € K, 

H(x,y) = H(x,z)H(z,y). 

Hence, there exists a signature function Sk on K [take H(-,zo) for a fixed 
zq in K] such that, for any x,y £ K, 

S K (x)S K (y)G(x,y) = \G(x,y)\. 

Denote by S n the function S[_ n n ] , then repeating the above argument and 
letting n tend to oo, we finally obtain the existence of S satisfying 

S(x)S(y)G(x,y) = \G(x,y)\. 

If S is not identically equal to 1, then there exists a point x at which S(x) = 1 
and liminfj^z S(y) = — 1. By continuity, this means that G(x,x) has to be 
equal to [because G(x,y n ) < for a sequence (y n )n>o that converges to x]. 
This is excluded since G is positive definite. Consequently, G is positive. 
We define the measure m on R by 



Note that the measure m is finite and / y/ G(y, y)m(dy) < oo. Making use of 
the covariance inequality (G(x,y) < \/G(x,x)G(y,y) for x,y in R), and the 
dominated convergence theorem, we see that the function 

(4) x(x) = / G(x,y)m{dy) 



is continuous. 

We now consider the fixed compact set K = [a, b] with a < b. For any 
integer n, we keep the definition of d n introduced at the beginning of the 
proof. We set 

I n = (k E N: -n2 n + 1 < k < n2 n - 1 and A G d n (if)|. 
Define G n on the set {( Jr, 4r) : k, i € I n } by 

Since the process rf is infinitely divisible and the positive measure m has 
a support equal to R, G" 1 is an M-matrix (G n need not be symmetric). 
Hence, we can write 

G n — c n I B n , 

where B n >0 and c n is strictly greater than the absolute value of any eigen- 
value of B„ . 
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Define 

Xn 



and let -D n be the diagonal matrix diag(x„( Jr), k G J n ). Then, -D n e n = G n e n , 
where e n is equal to (1, 1, . . . , 1). Consequently, for any k € I n , 



Setting T n = ^-D„ 1 B n D n . we have 

(5) D~ 1 G~ 1 D n = c n {I-T n ) 

and for every fe, X^£/ n T n (k,£) < 1. Consequently, the matrix T n is the tran- 
sition matrix of a transient Markov chain. 
Rewriting (5), it follows that 

(6) G n = D n O n (D n )-\ 

with O n = f (/-Tn)" 1 . 

Let ^4 be a square matrix of size |/ n | defined by A = (A(Jjfr, ^))k,lel n ■ We 
associate with A an operator on the set of the continuous functions with 
compact support. We denote this operator A and define it as follows. Let / 
be a continuous function on K with a compact support, then the function 
A/ is given by 

A/(s) = £ A[d n {x)^f[^j V.ER, 

with the convention that A(^c, ^r) = when k is outside of I n . 

That way we associate with D n (resp. D~ l , G n , O n ) the operator D n 
(resp. D~ x , G n , O n ). Note that we have, for every function / and every x 
in R, 

D„/(x) = Xn(^„(a;))/(d„(a;)) 
(D n )-V(x) = { Xn {d n {x))r l f{d n {x)) 

"^'^(wsj'yt m( ' i!/, • 

By (5), we know that the sequence (O n ) n >o satisfies 

O n = D n 1 G n D n . 

Moreover, for each n, O n is the Green operator of a finite state space Markov 
process. 
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Let O be the operator defined on Ck (the continuous functions with sup- 
port in K) by 

Of{x)=D- 1 gDf{x), x£K, 

where 

Df(x) =x K (x)f{x), 
D- 1 f(x) = ( XK (x)r 1 f(x) 

and 

Qf{x) = \ G(x,y)f(y)m(dy), 
Jk 

with 

X K {x) = \ G(x,y)m(dy). 

JK 

Note that, by the continuity of G and the fact that G(x, x) > and m has K 
as its support, x K ^ s continuous, strictly positive on K and, thus, bounded 
below by a strictly positive constant. 

Lemma 4.2. For every function f € Ck, the sequence (O n /) n >o con- 
verges to Of uniformly on K . 

Proof. For any continuous function / on Ck and any x 6 K, we have 

°nfi x ) = — ,] , u / G(d n (x),d n (y))xn{dn{y))f{dn{y))m{dy) 
Xn{d n {x)) Jk 

and 

Of(x) = — X—r I G(x,y)x K (y)f{y)m( d y)- 

X K \ x ) Jk 

Note that 

Xn(d n (x))= G(d n (x),d n (y))m(dy). 
Jk 

Since G is uniformly continuous on compacts and m is finite, the sequence 
(Xn(d n (x))) n£ ^ converges uniformly on K to x K ( x )- Hence, the sequence 
G(d n (x) , d n (y))xn (d n (y))f (d n (y) ) converges uniformly on K x K to G(x, y) x 

x K {y)f(y). □ 

We would like to show that O is a Green operator. For this, we shall use 
the following lemma which, in essence, is due to Hunt (see [12], Chapter X, 
page 255, or [15]). 
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Lemma C. Let V be a kernel on a measurable space (E,£) such that V 
satisfies the complete maximum principal and the function VI is bounded. 
There exists then a sub-Markovian resolvent (V p ) such that Vq = V . 

As Meyer has noted in [12], page 253, if we assume that V is continuous, 
then it is sufficient to verify the complete maximum principle for continuous 
functions with compact support only. More precisely, we have to verify that, 
for any a > 0, any two positive continuous functions with compact support, 
(/, h), if for all x in {h > 0} 

a + Vf(x)>Vh(x), 
then the inequality remains valid for all x in E. 

Lemma 4.3. The kernel O satisfies the complete maximum principle on 

K. 

Proof. First note that O maps continuous functions with support in 
K to continuous functions and is therefore a continuous kernel on K. Let / 
and h be in and a > and suppose that 

(7) a + Of{x)>Oh(x) for all x G {h > 0}. 

Recall that {h > 0} is contained in K. We need to show that (7) is satisfied 
by all x G K. Suppose that this is not true. Then there exists a constant 
b > and a point xq in K so that 

a + Of(x ) <Oh(x ) -b. 

Let e > and N be such that, for any re > N , 

|O n /(a?)-0/(x)|<e VxG/f 

and 

\O n h(x) -Oh{x)\ <e VxGif. 

Such N exists by Lemma 4.2 and the fact that /, h G Cj£. 
By (7), for every n > N and x in {h > 0}, 

(8) a + 2e + O n f(x) > O n h(x). 

Recall now that, for each n, O n is a potential of a Markov process on I n . 
It follows from [16] that O n satisfies the complete maximum principle on I n 
and, by its definition, O n satisfies the complete maximum principle on K. 
Consequently, (8) is valid for all x in K. In particular, we have 

a + 2e + O n f(x ) > O n h{x ). 
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On the other hand, we have 

a - 2e + O n / (so) < O n h(x ) - b. 
Choosing e < 6/4 leads to the desired contradiction. □ 

In order to simplify the rest of the proof of Theorem 3.4, we use now, 
instead of O, the operator Vk defined by 

V K f{x) = —^( G(x,y)f(y)m(dy). 



Xk 0) JK 

First we note that = 1. Further, since Vicf = 0(f /Xk)> Vk satisfies 
the complete maximum principle on K, by Lemma 4.3. 

Let K n be the compact set [— n, n] for n G N, and denote by V K the 
corresponding operator as defined above for K. Further, arguing as for x K > 
one can show that x defined in (4) is strictly positive. 

Lemma 4.4. Let V be the operator defined on bounded Borel functions 

by 

(9) y/( x ) = _L f G (x,y)f(y)m(dy). 

X{x) Jr 

Then there exists a sub-Markovian resolvent (V p ) such that Vq = V . 

Proof. Note that VI = 1, and that V maps continuous functions with 
compact support to continuous functions on R. Let / and h be two positive 
continuous functions with compact supports and a > and suppose that 

(10) a + Vf{x)>Vh{x) for all x G {h > 0}; 

we need to show that this is satisfied for all x in R. 

We denote by H a compact set that contains both the compact supports 
of / and of h. There exists no such that, for n > no, H is contained in K n . 
Hence, (10) can be written as 

a jdx) + y K j [x) > V Kn h{x) for all x G {h > 0}. 

Xli n \ X ) 

We have shown, when defining m, that x is a continuous function. Moreover, 
for any n G N, the function x Kn ^ s continuous and the sequence (x Kn ) * s m " 
creasing and converges simply to x- Consequently, by Dini's theorem, (x K ) 
converges uniformly to x an d x( x ) is strictly positive. Since the sequence 
(X K )n>no is bounded below by a strictly positive constant on H, x/x Kn 
converges to 1 uniformly on H. Hence, for every e > 0, there exists N such 
that, for every n> N, 

a(l + e) + V K J(x) > V K h{x) for all x G {h > 0}. 
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Since the operator V K satisfies the complete maximum principle on K n , the 
above inequality is still true for x in K n . We now multiply each side of the 
inequality by x K ( x )/xi x ) to obtain 

a(l+e)x Kn (x)/x(x) +Vf(x)>Vh(x) for all x G K n . 

Since x K J x )/x( x ) < 1, we finally get 

a(l + e) + Vf(x) > Vh{x) for all x G K n , 

and letting n tend to oo, 

a(l + e) +Vf(x)>Vh(x) for all s G R. 

Since this is true for any e > 0, (10) is established for all i in M. □ 

The sub-Markovian resolvent (T^,) allows one to construct a semi-group 
(Pt) with (V p ) as its resolvent. It will be a Ray process if we restrict our state 
space to a compact set, and, in general, we may need to apply a Ray Knight 
compactification in order to obtain a good semigroup. With this semigroup 
one can construct a transient Markov process on E = R with 0-potential 
equal to V. On the other hand, (9) implies that m is a reference measure for 
the Markov process with potential V and that the potential density h(x,y) 
of V with respect to m is equal to 

X( x ) 

Setting then fJ,(dy) = x(y) m (dy), we see that the operator V admits the 
symmetric densities (-^^G(x,y)^^,x,y G R) with respect to [i. □ 

Remark 4.5. Assume that G is such that, for any p, the inverse of the 
covariance matrix of (j]x\ > Vx2 > • • • i Vx p 

) is a diagonally dominant M-matrix. 
One can then follow the steps of the proof of Theorem 3.4, without the need 
to define the matrices D n (we check similarly, as in Remark 4.1, that the 
matrix T n = —B n is a transition matrix) and conclude that there exists a 
symmetric Markov process with potential densities equal to (G(x, y), (x, y) G 
R 2 ). This, together with Theorem 3.6(i), leads to Theorem 3.6(h). 

Finally, every centered Gaussian process with infinitely divisible square is 
equal to a deterministic function times a Gaussian process that is associated 
with a Markov process. The isomorphism theorem of Dynkin can hence be 
easily adapted to incorporate the deterministic functions, so that it can be 
used in studying any Gaussian processes with infinitely divisible square. 



INFINITELY DIVISIBLE SQUARED GAUSSIAN 
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5. The case of the Brownian sheet. A Brownian sheet is a centered Gaus- 
sian process (W x , s ,% > 0,s > 0) with a covariance given by 

E(W x , s W yjt ) = {x Ay)(s At). 

Remember that, up to a multiplicative constant, {x A y; x, y £ R+) is the 
Green function of the linear Brownian motion killed at its first hitting time 
of 0. Hence, this covariance is the product of two Green functions and one 
may ask whether (W^ s ,x > 0,s > 0) is infinitely divisible. The answer is 
given in the following proposition. 

Proposition 5.1. (i) For every (xi, Si)i<i< 3 the vector (W% Si , 1 < 

i < 3) is infinitely divisible. 

(11) The process (W£ s ,x > 0, s > 0) is not infinitely divisible. 

Proof. Using Griffiths criterion (or Bapat's criterion), it is easy to 
check that a sufficient condition for the infinite divisibility of the square 
of a three-dimensional Gaussian vector, indexed by {1,2,3} and with a co- 
variance g, is 

(11) g(i,j)g(k,k)>g(i,k)g(j,k), 

for any choice of k in {1, 2, 3}. 

Since any Green function satisfies (11), so does the covariance function of 
the Brownian sheet, as a product of two Green functions. 

By Bapat's criterion, we know that a covariance matrix G, such that 
G > 0, corresponds to an infinitely divisible square Gaussian vector if G^ 1 
is an Af-matrix. We choose (xj,Sj)i<j<4 such that < xi < X3 < X2 < X4 
and < S4 < si < S3 < S2- Let G be the matrix {{xi A xj)(si A Sj))i<«j<4- 
We compute the coefficient Gi\: 

G^2 = ^1^3S4(^2 - X?,){SZ ~ Sl)s 4 > 0. 

Hence, G -1 is not an M-matrix and the corresponding Gaussian vector does 
not have an infinitely divisible square. □ 

Acknowledgments. We would like to thank Gennady Samorodnitsky and 
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